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A region of the bovine coronavirus (BCV) genome that functions as a packaging signal has been cloned. The 
291-nucleotide clone shares 72% homology with the region of mouse hepatitis coronavirus (MHV) gene lb that 
contains the packaging signal. RNA transcripts were packaged into both BCV and MHV virions when the 
cloned region was appended to a noncoronavirus RNA. This is the first identification of a BCV packaging 
signal. The data demonstrate that the BCV genome contains a sequence that is conserved at both the sequence 
and functional levels, thus broadening our insight into coronavirus packaging. 


The coronavirus genome is a single-stranded, positive-sense, 
27- to 31-kb RNA molecule, the largest among all RNA vi¬ 
ruses. The genome is encapsidated by multiple copies of the 
nucleocapsid (N) protein and is packaged as a helical nucleo- 
capsid in the mature enveloped virion (for reviews, see refer¬ 
ences 14 and 25). During infection, coronaviruses synthesize 
six or seven subgenomic mRNAs which share 5' and 3' ends 
that are identical to the genomic RNA (for a review, see 
reference 29). 

How coronaviruses recognize the viral RNA(s) to ensure 
specific encapsidation and packaging is only beginning to be 
understood. Thus far, a packaging signal has been identified 
only for mouse hepatitis coronavirus (MHV). By using MHV 
defective interfering (DI) RNAs, a 69-nucleotide (nt) packag¬ 
ing signal that maps approximately 20 kb from the 5' end of the 
genome within the gene lb open reading frame was identified 
(Fig. 1) (8, 16, 28). Inclusion of the packaging signal in a 
non-MHV RNA is sufficient to allow the RNA to be packaged 
into MHV virions (30). Mutagenesis of the predicted bulged 
stem-loop structure of the MHV packaging signal disrupts the 
ability of the sequence to function as a packaging signal (8). It 
was also shown that an MHV strain A59 subgenomic DI RNA 
is packaged into virions when it contains the DI packaging 
signal, but the subgenomic RNA was reported to be packaged 
less efficiently than its DI genomic RNA (4). 

Sequences that are required for RNAs to be packaged by 
bovine coronavirus (BCV) have not been determined. BCV 
belongs to the same antigenic subgroup as MHV, and data 
based on the 3' one-third of the genome that has been cloned 
and sequenced indicate that the two viruses are closely related 
(1, 2, 13, 15). However, multiple pieces of data suggest that 
BCV is different from MHV with regard to RNA packaging. 
The first is that BCV packages subgenomic RNAs in addition 
to genomic RNA, whereas MHV packages very little, if any, of 
its subgenomics. BCV N and M mRNAs were shown to be 
more abundant in virions, on a molar basis, than was genome 
(11). The second difference is that the BCV defective RNA 
Drep (Fig. 1), which is replicated and packaged by BCV, does 
not contain any gene lb sequence (6). 

To gain further insight into the requirements for coronavirus 
packaging, we asked whether the BCV genome contains se¬ 
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quence information equivalent to the MHV packaging signal. 
Reverse transcription followed by PCR (RT-PCR) was carried 
out by using BCV virion RNA as the initial template. RT-PCR 
products were generated with primers 5' EcoRI (5'-GCA GA 
ATTCAATGGCGTAGTGG-3') and 3' BamHI (5'-GGCCT 
CGGATCCATACTTCTGAAT-3'). Primers were designed 
based on published sequence information for MHV (5) and 
unpublished sequence information from D. Yoo, University of 
Guelph, Guelph, Ontario, Canada. EcoRI and BamHI restriction 
sites (underlined) were included for subcloning. Multiple RT- 
PCRs were performed, and the products were pooled prior to 
isolation of the 291-nt DNA product, digestion with EcoRI and 
BamHI, and ligation into EcoRI/BamHI-cut pGEM3Zf(+). 
Four clones (designated pBCVpkgl to pBCVpkg4) were se¬ 
quenced by using Sequenase (USB). All four clones had the 
same sequence. The sequence revealed that the BCV genome 
contains a region that is homologous to the part of MHV gene 
lb where the packaging signal maps (Fig. 2). Alignment of the 
BCV and MHV sequences, using the Gap program in the 
Genetics Computer Group software package (27), indicated 
that the regions have 72% homology overall. The 69 nt of the 
BCV cDNA that aligned with the MHV packaging signal is 
74% conserved (underlined region in Fig. 2). 

After determining that a sequence exists in the BCV genome 
that shares homology with the MHV packaging signal, we 
analyzed the RNA sequence to determine its predicted sec¬ 
ondary structure. We were interested in determining the pre¬ 
dicted structure of the putative BCV packaging signal, since it 
was already demonstrated that the structure of the MHV pack¬ 
aging signal is important for its function (8). We used mFOLD 
3.0, the most current version, with parameters set at 37°C, 1 M 
NaCl, and no divalent cations for our analyses (17, 32). Anal¬ 
ysis was initially done in the context of a 190-nt region (nt 42 
to 232 in Fig. 2) that corresponds to the region of MHV that 
was previously examined (8). Two structures with overall free 
energies of —56.2 and —51.7 kcal/mol were predicted for the 
BCV 190-nt RNA (Fig. 3). We also examined the 69-nt region 
in the context of the entire 290-nt cloned region. Identical or 
similar structures were predicted for the sequence in the con¬ 
text of the larger RNA. Fifteen of the 23 structures predicted 
for BCV were identical or similar to the one shown for BCV in 
Fig. 3A. 

Parallel analysis of the 190-nt homologous region from 
MHV predicted only one structure (data not shown) that was 
similar to the predicted BCV bulged stem-loop with the lower 
free energy (Fig. 3A). However, the predicted structure was 
different from the structure that was reported previously for 
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FIG. 1. Comparison of the BCV and MHV strain A59 genomes and repre¬ 
sentative defective genomes. Drep is a cloned BCV defective RNA (6), and 
MIDI-C is a cloned MHV strain A59 DI RNA (7). The various parts of the 
genome that are included in each defective genome are indicated by dashed lines 
extending from the schematic of the parental genome. Drep includes a 30-nt 
reporter sequence (T) derived from the TGEV N gene. MIDI-C includes a 
known packaging signal (P). The approximate location of the packaging signal in 
the MHV genome (P) and the assumed position of the 291-nt cloned region of 
BCV (?) are indicated. Schematic representations of the 291-nt BCV RT-PCR 
and 312-nt MHV PCR products included in plasmids pBCVpkg and pMHVlb, 
respectively, are shown. NCR, noncoding region. 


the MHV packaging signal (8). This was surprising, since the 
reported structure was at least partially confirmed by muta¬ 
tions in base-paired regions that either disrupted or restored 
the packaging signal function (8). What accounts for this dif¬ 
ference is not clear. It is possible that changes in parameters 
used to predict RNA secondary structures account for the 
difference between the previously reported structure and our 
observations. It should be stressed that the secondary struc¬ 
tures shown here are only speculative in the absence of exper¬ 
imental data. Further experimental analysis is required to de¬ 
termine the validity of the predicted BCV structure(s). 
Nonetheless, the conservation of the BCV and MHV struc¬ 
tures predicted by our analysis justified further experiments to 
determine if the newly cloned BCV sequence contains a func¬ 
tional packaging signal. The results also suggested that the 
putative BCV packaging signal has the potential to fold into 
different structures. 

To determine directly whether the region we cloned from 
BCV could function as a packaging signal, Bpkg-CAT.R was 
constructed by appending the chloramphenicol acetyltrans- 


A 

C A 


NOTES 581 


AU 

GC 

A 

A G 

a ua g 

jX 

C A 
G A 
AU 
CG 

gUA a 
G r C 

A GC 

UG 

C u a 

G au a 

CG 

AU 

GC 

a au 

88 

gu a 

A 

AU 

CG 

UA Hr, 


A 

A 


-56.2 kcal 


U 

c u c u 

UG 

GC 

UA 

A A 
GC 
GC 

AU c 

C g 

H A C 
GU 
G A 

an 

AU 

aB 

UA 
GC 
5 ’ 3 ’ 


A 


’c g a 


U C°Cr 

G r 

U A A 

G 

A » C 

C ,rA 

U G C A 
_^ A 

A ! 
GC A 


G A 

AGU c G - 


GC 

Gc 

UG 


X 


CuAUCg 

X A A A 
\ AU 
CG 
AU 
GC 
CG 

u UG 

a au 

88 

gu a 

au a 

CG 

/X 

U r U 
AU *- 

68 

k 

A A 
GC 

X6 C 

c c 
c C 
UA 
GC 
GU. 

G, A 
AU 
GC 
AU 
AU 
AU 

ua 

c cc A 
UA 
G A 
UA 
CG 
5 ’ 3 ' 


-51.7 kcal 


FIG. 3. Predicted secondary RNA structures of the cloned BCV region. 
Predictions were made in the context of both the 190-nt (nt 42 to 232 in Fig. 2) 
and 291-nt (nt 1 to 291) RNAs. Only nt 59 to 200 (Fig. 2) are shown for 
simplicity. Arrows indicate the position of the 69-nt region that shares homology 
with the MHV packaging signal. 


ferase (CAT) gene at the 3' end of this sequence (Fig. 4). To 
accomplish this, the hepatitis delta virus ribozyme and T7 
terminator from plasmid v2.0 (21) were subcloned downstream 
of the CAT gene in pcDNA3.1/Zeo/CAT (Invitrogen) to yield 
pCAT.R. The CAT-ribozyme-terminator cassette was then re¬ 
leased from pCAT.R by cutting with Xhol and Sphl and sub¬ 
cloned into plasmid pBCVpkg. The CAT cassette was also 
cloned downstream of the MHV packaging signal in pMHVlb 
to generate Mlb-CAT.R (Fig. 4). Plasmid pMHVlb contains a 
312-nt fragment that was PCR amplified from MHV DI 


1 AAT G GCG TAG T GGT GGACAAG GT T GGAGACACAGAT TGTGTGTTTTATTT 
T G A TT T GTG AA GG 

51 TGCTGTGCGTAAAGAGGGTCAGGATGTCATCTTCAGCCAATTCGACAGCC 
C G C T GTACA GG 

101 T GAGAGT C AGC T CTAACCAGAG C C CACAAGG TAAT C T G GGGAGTAAT GAA 
TGA CCG CA T G C TG CGC 

151 CCCGGTAATGTCGGTGGTAATGATGCT CTGGCAACCTCCACTATCTTTAC 
GTG G C A C A A GCGTGG 

201 ACAAAGCCGTGTTATTAGCTCTTTTACATGTCGTACTGATATGGAAAAAG 
T A AT AT ATC C CC AT A G G 

251 AT T T TATAG C T T TAGAT CAAGAT GT G T T TAT T CAGAAG TAT 
G A G T C GCA A 
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FIG. 2. Sequence of the cloned BCV region and alignment with the homol¬ 
ogous region from the MHV strain A59 genome. The region of the BCV clone 
that is homologous to the MHV packaging signal is underlined. Only noniden¬ 
tical nucleotides are shown for MHV. 


FIG. 4. Schematic representation of constructs used to test the ability of the 
BCV clone to function as a packaging signal. All coronavirus sequences were 
placed under the control of the T7 promoter. A cassette containing the CAT 
gene, the hepatitis delta virus ribozyme (r), and the T7 terminator (T7t) was 
cloned at the 3' end of each coronavirus sequence. 
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MIDI-C (7) and subcloned into .EcoRI and Bam HI sites of 
pGEM3Zf(+). This construct was used as a control for pack¬ 
aging in MHV-infected cells and to determine if the BCV and 
MHV sequences are functionally equivalent. Two additional 
constructs, BN3'-CAT.R and MN3'-CAT.R, which had been 
made previously for other reasons were also used as controls in 
the experiments described below. These included an 830-nt 
region and a 385-nt region from the 3' ends of the BCV and 
MHV N genes, respectively, upstream of the CAT cassette 
(Fig. 4). 

To determine if the addition of the cloned BCV sequence 
was sufficient for packaging of the CAT gene, we performed 
the following experiments. BHK cells for the BCV experiments 
or 1701 cells for the MHV experiments were infected with 
vTF7-3 (9), the vaccinia virus recombinant that expresses T7 
RNA polymerase. After infection, cells were transfected with 
the different plasmids essentially as described previously (19) 
(Fig. 4). Plasmid pMHV4S, which contains the MHV4 spike 
(S) gene (10), was included in transfections in the BCV exper¬ 
iments since fusion of the monolayer by the expressed S pro¬ 
tein helped increase the number of dually infected and trans¬ 
fected cells. Four hours after transfection, cells were either 
mock infected or infected with BCV or MHV. At 24 h after 
infection with BCV or 12 h after infection with MHV, both 
intracellular RNA and the media were harvested. The media 
were clarified and treated with both DNase and RNase prior to 
being either pelleted through a 30% sucrose cushion or sedi¬ 
mented on a 20 to 60% sucrose gradient. After gradient frac¬ 
tionation, the fractions with densities ranging from 1.15 to 1.20 
g/cm 3 were pooled and virions were pelleted. RNAs were ex¬ 
tracted from pelleted virions. One-tenth of the total intracel¬ 
lular RNA (Fig. 5A and B, lanes 1 to 10) and all of the RNA 
from pelleted virions (Fig. 5A and B, lanes 11 to 20) were 
analyzed by Northern blotting using a 32 P-labeled CAT-spe¬ 
cific riboprobe. Densitometry was used to quantify bands after 
autoradiography of Northern blots. After autoradiography, the 
membranes were reprobed with N gene-specific probes to con¬ 
firm successful coronavirus infection and the absence of coro- 
navirus infection in the controls (data not shown). 

Northern blotting clearly demonstrated that RNA tran¬ 
scribed from our cloned BCV sequence was sufficient to target 
CAT for packaging into extracellular BCV virions (Fig. 5A, 
lane 16). The two bands detected by the CAT-specific probe 
likely resulted from incomplete ribozyme activity. Interest¬ 
ingly, the data also showed that the BCV RNA is functionally 
homologous to the MHV RNA since it was also sufficient for 
packaging of a non-MHV RNA into MHV virions (Fig. 5B, 
lane 18). Consistent with previously published data, the MHV 
packaging signal was functional in MHV-infected cells (Fig. 
5B, lane 16) (4, 30). In addition, the results provided new data 
indicating that the MHV signal is also functional in BCV- 
infected cells (Fig. 5A, lane 18). In the case of extracellular 
BCV virions, small amounts of both the BN3'-CAT and CAT 
transcripts were detected in virions in some experiments (Fig. 
5B, lanes 17 and 19). Only very small amounts of CAT and 
MN3'-CAT were detected in extracellular MHV virions in 
some experiments (Fig. 5B, lanes 17 and 19). 

The ratio of each packaged RNA relative to the amount of 
its intracellular RNA was determined by densitometric scan¬ 
ning of RNA bands on Northern blots from multiple experi¬ 
ments like those shown in Fig. 5. The average of multiple 
experiments indicated that the Bpkg-CAT and Mlb-CAT 
RNAs were packaged into BCV virions about five- and sixfold, 
respectively, better than the CAT or BN3'-CAT transcripts. 
Both RNAs were packaged about sixfold better in MHV strain 
A59-infected cells. 
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FIG. 5. Northern blot analysis of expressed RNA transcripts packaged by 
BCV (A) and MHV strain A59 (B). Plasmid DNAs were transfected into vTF7- 
3-infected cells (lanes 2 to 9 and 12 to 19). Both intracellular RNAs (lanes 1 to 
10) and extracellular virion RNAs (lanes 11 to 20) were analyzed by Northern 
blotting after cells were mock infected (lanes 2 to 5 and 12 to 15) or infected 
(lanes 6 to 9 and 16 to 19) with BCV (A) or MHV strain A59 (B). Extracellular 
media were treated with both DNase and RNase prior to isolation of virions. 
RNAs were analyzed by using a CAT-specific riboprobe. Arrows indicate the 
position of the packaged Bpkg-CAT and Mlb-CAT RNAs. 
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Taken together, the results of this study demonstrated for 
the first time that the BCV genome contains a sequence that 
functions as a packaging signal, at least when appended to a 
noncoronavirus RNA. The BCV signal exhibits 74% homology 
with the MHV packaging signal. This conservation is sufficient 
to allow the two sequences to be reciprocally functional. It 
remains to be determined whether the sequence we have iden¬ 
tified, like that of MHV, is the functional packaging signal in 
the context of the wild-type BCV genome. For both viruses, 
this must await the construction of an infectious clone. 

The MHV packaging signal is necessary for packaging of 
several DI RNAs (8, 28). Since the BCV defective genome 
Drep does not contain any of the genomic sequence we have 
identified here, it would be informative to determine if the 
BCV signal increases the packaging efficiency of Drep when 
incorporated into the defective RNA. Thus far, our attempts to 
assay the effects of incorporation of the sequence into Drep 
have not been successful. We have introduced the sequence in 
frame into multiple sites within the defective genome, and in 
all cases replication of the defective genome was completely 
compromised (R. Cologna and B. G. Hogue, unpublished data). 

Since BCV and MHV contain highly homologous packaging 
signals, why does BCV not discriminate among the different 
viral RNAs whereas MHV is more selective in packaging? 
Since BCV packages subgenomics and the Drep defective ge¬ 
nome, none of which contain the identified BCV packaging 
signal, other factors or sequences may function as alternative 
signals for these RNAs. Alternatively, subtle differences be¬ 
tween the nucleocapsid proteins of the two viruses may ac¬ 
count for the packaging difference. In the case of alphaviruses, 
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the ability to package genomic RNA over the subgenomic 
RNAs was lost when a small region of the capsid protein was 
deleted (20). The nucleocapsid interacts with genomic and 
subgenomic RNAs in both MHV-infected cells (3, 26) and 
BCV-infected cells (R. Cologna and B. G. Hogue, unpublished 
data). The BCV N-subgenomic RNA interactions may be more 
stable than the comparable MHV complexes, thus allowing 
BCV subgenomics to be readily incorporated into assembling 
virions. Another possibility is that replication and assembly 
complexes are compartmentalized in MHV-infected cells, 
whereas these complexes may not be separated in BCV-in¬ 
fected cells. This might allow the abundant subgenomic RNAs 
to be nonselectively included during assembly. Our data sug¬ 
gest that BCV is less selective than MHV in packaging since a 
small amount of the CAT RNA was consistently detected in 
BCV virions. Finally, differences in the translation efficiencies 
of subgenomic RNAs from the two viruses might account for 
which RNAs are packaged. For example, if MHV mRNAs are 
translated more efficiently than BCV mRNAs, more of the 
subgenomics in MHV-infected cells may be associated with 
polyribosomes, resulting in their being unavailable for packaging. 

Other coronavirus DI RNAs are currently being studied, but 
no packaging signal for them has been definitively identified. 
Defective RNAs have been examined for the avian infectious 
bronchitis virus, an antigenic group III virus. The data indicate 
that information at the 5' end of the lb gene may be involved 
in the packaging of these defective RNAs (22, 23). DI RNAs of 
transmissible gastroenteritis coronavirus (TGEV), a member 
of antigenic group I, have also been isolated and studied (12, 
18). Recent analysis of one of the DI RNAs suggests that two 
regions play a role in packaging (12). These regions include 
about 1 kb from the 5' end of the genome which consists of 
sequences from both genes la and lb in the DI RNA and 
about 4.1 kb at the 3' end of gene lb. The exact role of the 
regions in packaging remains to be shown. Both TGEV and 
infectious bronchitis virus package subgenomic RNAs (24, 31). 
Clearly, further studies are necessary to comprehensively un¬ 
derstand the requirements for coronavirus packaging. Our 
data provide a step in that direction. 

Nucleotide sequence accession number. The packaging sig¬ 
nal sequence was deposited in GenBank under accession no. 
AF159126. 
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